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Abstract 

We give a nonperturbative derivation of the Bethe-Salpeter equation based 

on the Feynman-Schwinger path integral representation of the one-particle 

propagator in an external field. We apply the method to the quark-antiquark 

system in scalar QCD and obtain a confining BS equation assuming the Wilson 

area law in the straight line approximation. The result is strictly related to 

the relativistic flux tube model and to the qq semirelativistic potential. 
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I. INTRODUCTION 



Various attempts have been done in the literature to apply the Bethe-Salpeter equation 
to a study of the spectrum and the properties of the mesons. The hope was to obtain an 
unified and consistent description of the quark-antiquark bound states envolving light quarks 
as well as heavy ones. In all such attempts, at our knowledge, the choice of the kernel of 
the equation was purely conjectural and only made in such a way that the successful heavy 
quarks potential could be recovered in the non relativistic limit. While a kind of derivation 
of the qq potential from QCD can be given in the Wilson loop context, even if at the price 
of not completely proved assumptions [0- ||, no similar derivation seems to exist for the 
BS equation. 

In this paper we want to consider this problem for the simplified model of the scalar 
QCD. 

As it is well known a satisfactory semirelativistic potential for the heavy quark-antiquark 
system can be obtained from the assumption 

i In W = f(ln W) pcvt + aS min (1.1) 

where: W stands for the Wilson loop integral 

W = l(TrP expig{f dafAJ) ; (1.2) 

T denotes a closed loop made by a quark world line (Ti), an antiquark world line (r 2 ) 
followed in the reverse direction and two straight lines connecting the initial and the final 
points of the two world lines; (lnK / ) pert and S m { n are the perturbative evaluation of In W and 
the minimum area enclosed by T respectively; finally the expectation value in (1.2) stands 
for the functional integration on the gauge field alone. More sofisticate evaluations of i In W 
have also been attempted 0, ||], J7], but they shall not be considered here. 

In the derivation of the potential, S m m is further approximated by the surface spanned 
by the straight lines connecting equal time points on the quark and the antiquark worldlines, 
i.e. by the surface of equation 
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x° = t, x = szt(t) + (1- s)z 2 (t) , (1.3) 

t being the ordinary time, Zi(£) and z 2 (t) the quark and the antiquark positions at the time 
t, and s a parameter with < s < 1 . The result is 

S mi n = £ dtr£ds[l - (sz 1T + (1 - s)z 2T ) 2 ]5 = (1.4) 

where Zj T stands for (S hk - r ^r)^ and r(t) = z 1 (t) - z 2 (t). 

In fact it can be shown that (1.4) is correct up to the second order in the velocities 
H, |§ and so is perfectly appropriate for a derivation of the semirelativistic potential. For 
a full relativistic extention, Eq.(1.4) cannot be correct as it stands, due to the privileged 
role played by the time in it. However one can try to assume (1.4) in the center of mass 
frame and then obtain implicitly the corresponding equation in a general frame simply by 
Lorentz transformation. In this way the assumption would become equivalent to the so 
called relativistic flux tube model 0,0. 

In this paper we want to show that in scalar QCD (i.e. neglecting the spin of the quarks) a 
Bethe Salpeter (BS) equation which include confinement can be derived from (1.1) and (1.4) 
in the center of mass frame. In the momentum representation and after the factorization of 
the four-momentum conservation 5, the resulting kernel turns out in the form 

I(Pl,P2\ PliP'2) = W(P1,P2! p'i,P2) + £anf(pi,P2; P'liP'2) ( L5 ) 

{p'i +p' 2 = Pi +P2), where i pert is the usual perturbative kernel which can be written, at the 
lowest order in the strong coupling constant, 

W(Pi>P25Pi»P2) = Y 9 ( — 2 — ) ^Vfri -Pi)( — 2 — ^ ' ^ ^ 

while in the center of mass system J con f is given by 

T 1 I f\ f j3 ifk'-kVr r( Pi + Pi P% + P2\ /i ^ 

IcoAPi,P2;Pi,P 2 ) = j d J re 1 >J{r, — - — , — - — ) (1.7) 
(pi = — P2 = k, p[ = — P2 = k') with J expressed as an expansion in 
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J(r, qi, 92) = 4^ [?2o\/9io - Qt + 9?o\M) - Qt + 

^ <?10 + 920 



,2 „2 



H — (arcsm- + arcsm- -)J + . . . (1.8) 

Iqt| |?io| mol 



Notice that by a usual instantaneous approximation, from the kernel defined by ( |1.5|) -(|1.8| 
one can obtain the following hamiltonian 



H{r, q) = y mf + q 2 + \Jm 2 2 + q 2 + 

or 1 J 

9 /„__2 1 19 1 /.._2 1 I9V 



- : 1 ^2 + q f—j 2+ "°i 1 * 

2 + g 2 + */m| + ^ m i + 9 2 ?r ^2 + 9 



m i + 5 / — 2 



+ q z 
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'mf + qym^ + q\ f . q T q T x-, 

1 I arcsm— 7= + arcsm— = )>• + ... 

<?T yv ^m 2 + q 2 ^m 2 + q 2J) 

+ V pcrt (r,q) (1.9) 

with an appropriate ordering prescription. In Eq. ( |1.9|) q stands now for the momentum in 
the center of mass frame, q r = (f • q)f, = (5 hk — f h r k )q k , while V per t is the ordinary 
perturbative Salpeter potential. We stress that Eq.( |1.9| ) is identical to the hamiltonian for 
the already mentioned relativistic flux tube model and consequently is strictly connected to 
the semirelativistic potential as given in [Q, || when the spin dependent part is neglected. 
Notice however that the ordering in (1.9) corresponding to (1.7) is not identical to the Weyl 
prescription. 

The present paper has mainly a pedagogical purpose and we have not done any attempt to 
explicitly apply the kernel (1.5)— (1.8) to an evaluation of the spectrum and of the properties 
of the mesons. Notice however that very interesting results have been obtained in this 
direction by the relativistic flux tube model || and in a sense we may consider our paper 
also as providing a more fundamental justification to that model. 

As in potential theory the starting object in our derivation is the gauge invariant quark- 
antiquark propagator 

Gf (xi,x 2 ; 2/1,1/2) = (0 I (j)* 2 (x2)U(x 2 ,xi)(f) 1 (x 1 )(f)* 1 (y 1 )U(y 1 ,y2)(t)2(y2) I 0) = 

= ~Tr(U(x 2 ,x 1 )A$\x 1 ,y 1 ,A)U(y 1) y 2 )AP(y 2 ,X2,A)) , (1.10) 



where U(b,a) = P& a exp(i<7 j a dx^A^x)) is the path ordered gauge string (the integration 
path is over the straight line joining a to b), while and Ap denote the Feynman 
propagators for the two quarks in the external field A^. In contrast with the potential case, 
however, no semirelativistic expansion is used, but the propagators are treated exactly using 
the covariant Feynman- Schwinger path integral representation. 

Notice that in establishing the BS equation we have to neglect in (ilnW) peTt the con- 
tributions from the two extreme lines x\x 2 and y\y 2 and in S m i n the border contribution 
corresponding to x w ^ x 2 o and y w ^ y 2 o, as in the potential case. This is correct for x® — y® 
and x 2 — 2/2 large with respect to |xx — x 2 |, |yi — y2 1 , — x 2 and — y 2 . So, strictly, we 
obtain a BS equation but for a quantity G4 which coincides with Gf in the above limit. 
This is immaterial for what concerns bound states or asymptotic states. Naturally out of 
the limit situation, G4 is no longer gauge invariant as it should be talking of a BS equation. 

The significance of the Feynman-Schwinger representation in the framework of QCD has 
already been appreciated in JED] and particularly in [|]]. With respect to Ref. ||] we have 
however a different attitude on the role of the BS equation. 

The plan of the paper is the following one: in Sec. 2 we illustrate our method on the 
example of a one dimensional particle in a velocity dependent potential. In Sec. 3 we apply 
the method to the derivation of the Bethe-Salpeter equation for a system of two scalar 
particles interacting via a scalar field and discuss the various complications related to the 
perturbative kernel. In Sec. 4 we derive the BS equation for a quark- ant iquark system in 
scalar QCD under the assumption discussed above and obtain the kernel reported. Finally 
in Sec. 5 we derive the Salpeter potential from the BS equation. 

II. ONE DIMENSION POTENTIAL THEORY 

Let us consider the model made by a nonrelativistic particle in one dimension with the 
Hamiltonian 

2 

H=^- + U(x,p) = H + U (2.1) 
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and the corresponding Schrodinger propagators 

K(x,y,t) = (x\e- lHt \y) , K Q (x,y,t) = (x\e~ iHot \y). (2.2) 

From the operatorial identity 

e~ lHt = e~ lHot - i f di!e lH ^-^Ue- %m ' , (2.3) 

JO 

we obtain the equation 



o 



K(x,y,t)=K (x,y,t)-i dt' / d£ / d V K (x,$,t - t)(£\ U\q) K( V , y,t'). (2.4) 



which is somewhat analogous to the nonhomogeneous Bethe-Salpeter equation in the con- 
figuration space. 

We want to derive Eq.( [2.4|) by the path-integral formalism. 

Let us be for definiteness 

U = V(x) + (W(x)p 2 ) OTd , (2.5) 

where ( ) or( j stands for some ordering prescription. In terms of path-integral we can write in 
the phase space 

K(x,y,t) = [ X VzVpexp{i f dt'lp'z' - | V(z') - W(z')p' 2 ]} (2.6) 

J y Jo 2m 

with z' = z(t'), p' = p(t'), z' = ^jjp-- In Eq. (^6|) the functional "measures" are supposed to 
be defined by 

Vz = {——) *dz x ... d^Ar-i, Vp = (- ) 2 dpi . . . dpN-xdpN 

2me ZTim 

VzVp = (— ) N dpidzi . . . dp N _idz N _idp N (2.7) 
2n 

where e — 4p, the limit iV — > oo is understood and the end points x and y stand for the 
condition z = y, zn = x. As well known (see e.g. ||) the ordering prescription is concealed 
under the particular discretization adopted in the limit procedure implied in the definition 

of (Mi. 



Let us first consider the case W = 0. In this case there is no ordering problem and 
it is possible to calculate explicitly the integral in p in ( |2.6[ ) obtaining the path-integral 
representation in the configuration space 

K(x,y,t)= f X Vze l Io dt '( m ^- v ^\ (2.8) 
Jy 

Then using the identity 

e -*/>'vV) = 1 _, rdt'Viz^So^") (2 .9) 



o 

one obtains 

K(x,y,t) = K (x,y,t)-i f dt' (* VzVUy& M '^ M ti ^^ti (2.10) 

JO Jy 

which, taking into account that 

Vz... = Jd£J VzJ Vz... (2.11) 

(having identified £ = z(t')), can be rewritten in the form ( |2.4| ) with (£|C/|ry) = V(^)S(^ — rj). 

In the general case W(x) ^ it is convenient to work with the original path-integral 
representation in the phase space (|2.6|) and it is necessary to use discretized expressions 
explicitly. For Weyl ordering in Eq. (|2.5| ) the correct discretization is the mid-point one and 
we can write 



K(x, y, t) = 1 J dp N dz N -idp N -x . . . dz x dp x 



(27T) 

ex.piJ2{Pn{Zn - Zn-l) " ^ + v f*^ 1 ) + W(^^)^}. (2.12) 

Then it is convenient to introduce the Fourier Transform of K(x,y,t) 

K(k,q,t) = J dxj dye' tkx K{x,y,t)e iqy (2.13) 
and to use the discrete counterpart of Q2.9| ) 

exp (-ie) ) + W{ K] = 1 - is 1^{[V( ) + W{ )p 2 R ] 

n=l / Z R=l A / 

■ exp (-is) J2[V(^±^) + W(^±^zl)p% (2.14) 

r=l Z Z 



Replacing ( |2.14| ) in (|2.12|) and Fourier transforming we have 

N „ A 

(27T)" 



I f 

K(k, q, t) = K (k, q, t) - is J dz N dp N . . . dp x dz e~ lkZN exp i Pn(z n ~ Zn-i) 



R=l n=R 
R-l 



\V{^^) + W( Z -^±p^)p R ] exp* £ \p n (z n - ^ 

Z Z n=l 

v 1 r r , r . 1 



ikzisi 



= K (k, q, t)-ieJ2 J2^y{ J d PR+^ J d PR-i ( 2ir )N-R-i J dz N d PN ■ ■ ■ dp R+2 dz R+1 e 

N 

exp{2 Pn(Zn ~ Zn-l)}e lPR+lZR+1 } 

n=R+2 

I dz R d PR dz R-^~ l{PR+1 ~ m)ZR \ V C R + 2 ZR ~ 1 ) + W( ZR + 2 2?fl " 1 )Pfl]e i(Pfl - 1 - pfl)Zfi - 1 } 
1 /• R ~ 2 

( /0 x B _ 9 / ^-2 • • • dZ e~ %PR - lZR - 2 exp Z V {^(^n - Zn-l) 

H 27r ) ^ „=i 

_ £ [|L + V( z _iL+^zl ) + H /(^±^)^]}e^} (2.15) 
which in the continuous limit reads 

K(k, q, t) = K (k, q,t)-i f dt' J ' ^ J ^K (k, k', t - t')I w (k\ q')K{q\ q, f) (2.16) 



with 



I w (k,q) = ±-J dz"dpdz'e- ikz "e ipz "{V(^^) + W(^^)p 2 ] e - ipz ' e tqz ' = 

= V{k -q) + W(k - g)(^) 2 = (k\(V(x) + l{p, {p, W(x)}})\q). (2.17) 

Eq. fl2.16| ) is equivalent to Eq.(2.4) and the kernel reduces to the Fourier transform of the 
case W(x) = 0. 

Had we considered the symmetric ordering in Eq. ( |2.5|) we should have replaced in (|2.14|) 
and (Hi [V(^§^) + W(^f^)p 2 n } with [\{V{x n ) + V{x n ^)) + \{W{x n ) + W{x n ^))pl\ 
and the result would be 

I s (k,q) = — J dz 'dpdz e %kz e %pz [— — — — - H — — — p 2 ]e ipz e* 9 * = 

= V(k -q) + W(k - q)^^ = (k\(V(x) + \{p\ W(x)})\q) (2.18) 



8 



III. SPINLESS PARTICLES INTERACTING THROUGHT A SCALAR FIELD 



Let us consider two scalar "material" fields <p\ and 4>2 interacting through a third scalar 
field A with the coupling + gi(j)\A). Then, after integration over §\ and 02, the full 

one particle propagator can be written as 

G f{x-y) = (01^,(^)0,(^)10) = (iA«\x,y;A)) es ^^^^^^ (3-1) 

where Ap (x, y, A) is the propagator for the particle j in the external field A, Sq(A) is the 
free action for the field A and the determinantal factor M(A) comes from the integration of 
the fields <fij 

,2 



_ „ Md^ + m)- gj A) 



( ] 2 det(d^ + m 2 ) 

= 1 - \ £ { - 9i j d 4 xA(x)A^(0) - l -g) J d 4 xd 4 yA(x)A ( i\x - y)A(y)A?(y - x) + . . .} (3.2) 
The covariant Feynman-Schwinger representation for Ap reads 

A^(x, y- A) = -- J o dr J VzVpexp {i ^ dr'[-p^ + -p'^ - -m 2 + -gjA{z')]} 

j poo px rr I 

= ~\ I dr J Vzexp { - i jf dr'-[(i' 2 + m 2 ) - ^A(z')]}, (3.3) 

where the path integrals are understood to be extended over all paths z^ 1 = z^(t') connecting 
y with x expressed in terms of an arbitrary parameter r' with < r' < r. In Eq. (|3.3| ) 
stands for z(r'), p' for p(r'), i' for and the "functional measures" are assumed to be 
defined as 

Vz = (J-y N d 4 Zl . . . d 4 z N ^, Vp = (^) 2N d 4 Pl . . . d 4 p N _,d 4 p N 

2ms 2tt 

VzVp = (—) 4N d 4 Pl d 4 Zl . . . (Ppn-kPzn-kPpn. (3.4) 

Replacing Eq. (|3T3"|) in ( p.l| ) we obtain 

Gf{x-y) = i /""dr f X Vzexp{- 1 - f dr\z' 2 + m 2 )}(exp *§> f dr' A{z')) (3.5) 
2 Jo Jy 2 Jo 2 Jo 

where, suppressing the tadpole term in (3.2), 



( exp % 1L f T dr'A{z'))= exp ^ f dr' f dT>'[D F {z' - z") + 
2 Jo 4 Jo Jo 

+ E f d 4 V D F (z' - 0(A«(e - r ? )) 2 D F (r / - z") + . . .] (3.6) 

Now let us consider the two particle propagator 
G A (x 1 ,x 2 ;yi,y2) = (O|T0i(xi)0 2 (x 2 )0i(2/i)0 2 (2/2)|O) = (G i 2 1 \xi,yi;A)G i 2 2) (x2,y 2 ;A)) = 

[ roc roc rx\ rX2 

= (-) 2 / dn / dr 2 / Vz\ \ Vz 2 x 

Z JO JO Jyi Jy 2 

exp^{/ Tl dr[(z' l 2 + m\) + H dr' 2 (z' 2 2 + m 2 )}(exp % -{ 9l H dr[A(z[) + g 2 H dr' 2 A(z' 2 )}) . (3.7) 
z Jo Jo z Jo Jo 



In this case 



(exp % -{ gi f 1 dr{A{z[) + g 2 H dr' 2 A{z' 2 )}) = 
2 Jo Jo 

E l 4 P dT 'i f* d ^ D ^ - 4) + E I / *Z / <*VM4 - o x 

.7=1,2 4 i=l,2 

(A?(£ - ??)) 2 ^ - z'l) + ...] + ^ T ^ P ^[D F (^ - 4) + 

4 jo jo 

+ E f / A / A^f(4 - 0(A?(£ - r / )) 2 D F (r / - + . . . ] (3.8) 

If we replace M(A) by 1 in ( |3.1|) and ( |3.7|) , i.e., if we retain only the lowest order terms in 
the exponent in ( |3.8|) (quenched approximation), we have no "material" fields loops in the 
evaluation of the gauge field average ( p.6|) and ( p.8|) and then we can write exactly 



1 9 f°° f°° f x l f x 2 

G 4 (xi,x 2 ;y u y 2 ) = (-) / dr x dr 2 Vz x / Vz 2 

I Jo Jo Jyx Jy2 

exp Z![ r dr^ 2 + m 2 ) - | / Tl ^ f* <D F (4 - 4')] x 

Z JO Z JO JO 

exp Zl[ P dr^(4 2 + m 2 )} - | P ^ P ^'D F (4 - 4')] x 

Z JO Z JO JO 

exp ^ P dr( P dr'DM - z'). (3.9) 
4 Jo jo 

Proceding as in Sec. 2, using the identity 

exp — — / dr x \ dr 2 D F (z 1 - z 2 ) = 
A Jo Jo 

= 1 + ^ P dr[ P dr' 2 D F {z[ - z' 2 ) exp[^ P dr'( P d^D F {z'{ - 4)] (3.10) 
4 Jo jo 4 Jo Jii 
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( analogous to Q2.9Q ) and Eq. , we obtain after some manipulations 



G A {x u x 2 \ y 2 ) = G 2 (x 1 - yi)G 2 (x 2 - y 2 ) + 



(-) 2 / dr[ \ dn / dr' 2 / dr 2 \ d A z[ / d A z' 2 / Vz x \ Vz 2 / Vz x \ Vz 2 D F (z[ - z' 2 ) 

2 JO Jt[ JO Jt' 2 J J Jz[ Jz' 2 Jy! Jy 2 

exp l -{ jT dr'l\-^ - m\ + | fj dr'l' D F {z'{ - z'{')] + jT dr''\-Q - m\ + | jf* < 

D F (4 - 4")]} x exp l -{ f 1 dr'{[-z[ 2 - m\ + f d7?' Z>p(s? - 4")] + f * < H? ~ m 2 + 

Z Jo I Jo Jo 

+| f*' AK4' - 4")] + f * <^(4' - 4')} x 

1 Jo 1 Jo Jo 



exp*{f / n dr'l f' 1 dT?D F (4 - z'{') + f f" dr? P d^'D F (z' 2 ' - 4") + 
4 Jt[ Jo 4 Jr^ Jo 

f' 1 dr'l r 'dT%D F (z'{-4)}. (3.11) 

JO Jt' 



+ 9i92 f 1 j.n 



Then, if we denote by L\ the last exponential in this equation and replace it by 1, we obtain 
immediately the Bethe-Salpeter equation 

G^(x 1 ,x 2 ,y 1 ,y 2 ) = G 2 {x 1 - yi)G 2 {x 2 - y 2 ) + 

-i j d A z 1 J d A z 2 J d A z[ J d A z' 2 G 2 {x l - z 1 )G 2 (x 2 - z 2 )I{z u z 2 ; z[, z 2 )G A (z[, z' 2 ,y 1 ,y 2 ) (3.12) 
with the ladder approximation kernel 

I(zt, z 2 , 4, 4) = ~9i92Df(z[ - z' 2 )S A (z 1 - z[)S A (z 2 - z 2 ). (3.13) 
On the contrary if we introduce in (|3.11| ) the entire expansion of L\ in the second line of the 



a 

equation beside D F (z[ — z 2 ) we obtain additional terms of the type =f f , dr-' JZ* dr"' D F (z" — 
z'")D F (z[ — 4), 9\92 It[ dr" Jq 2 D f (z" — z 2 )D F (z[ — z 2 ) etc.. As a consequence (see App.A 
for details) we reobtain Eq. fl3.12|) but with the kernel 

I(zi, z 2 , z[, z' 2 ) = 

= -g l92 D F (z[ - z' 2 )5 A (z[ - z 1 )S\z' 2 - z 2 ) + ig\g 2 J d%D F ( Zl - z[)G 2 ( Zl - £i)G 2 (6 - 4) 
Df{£i - z 2 )5 A (z 2 - 4) + i 9l glj d£ 2 5 A ( Zl - z[)D F ( Zl - &)G 2 (z 2 - 6)G 2 (6 - 4) 
D F (z 2 - 4) + %g\g\D F {z x - z' 2 )G 2 { Zl - z[)D F (z 2 - z[)G 2 {z 2 - 4) + • • • (3.14) 
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or, 



graphically: 




Finally to go beside the quenched approximation and to take into account additional 
terms in Eq. (|3.1| ), Q3.8|) amounts to insert 4>i4>i and 4> 2 (f) 2 loops in all possible way inside the 
graph. 

Had we used Eq. ( |3~5|) in the phase space and a method analogous to that employed start- 
ing from Eq.(2.12) we would have obtained the Bethe-Salpeter equation in the momentum 
space 



G f 4 (pi,P2;pi,P 2 ) = G { 2 \pup'i)G 2 2 \p2,p' 2 ) 
d 4 k[ d A k' 2 d A h d 4 k 2 ~ 



(2vr) 4 (2nY (2tt) 4 (2tt) 



G 2 ( P i, h)G 2 (p 2 , k 2 )I(h, h; K, k' 2 )G 4 (k[, k' 2 -p' llP ' 2 ). (3.15) 



i 



with 



(2tiY (2tt) 4 

Agig2 D{^±A _ Z _1±A) jU-WjU-W (3.16) 

where for defmiteness we have assumed the mid-point prescription in the discretized form 
of (3.9) even if immaterial in this case. Then, introducing the total momentum P = P\ + p 2 
and the relative momentum q = — ^ — pi r r 1 — r> 2 , defining 

G 2 (p,p') = {2nY5\p-p')G 2 {p) G^p^p'^) = (2vr) 4 5 4 (P - P>)G±(q, q'; P) 

I(pi,P2; P [,P^ = (27r) 4 5 4 (P-P')i(q,q , ;P) (3.17) 

and factorizing the total momentum conservation delta, we can write in conclusion 

G 4 {q-q\P) = {2*)H\q-q')G%\q)Gf{-q) -i J -^G { 2 1 \q)G 2 2 \-q)i(q,k,P)G,(k,q',P) 

(3.18) 
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and in the ladder approximation we would have obtained 

i(q,q'; P) = - 9l g 2 D F (q - q') 



(3.19) 



IV. SCALAR QCD 



Let us come to the scalar QCD characterized by the lagrangian 

L = {D p 4>yD?4> - m 2 ct>*4> - ^Tt V" + ^gf (4.1) 

where D p = d p + ieA p and Lqf is the gauge fixing term. In this case an equation analogous 
to (13.31) can be obtained 



A$\x,y;A ll ) = --P xu J o dr J Vzex V {-i jf dr'[-(z' 2 + m) ) - z^A^z')}} 
A^(y,x;A p ) = --P yx dr Vzexp{-i dr'[-(z' 2 + m 2 ) + z* A p (z'))} . (4.2) 

Z JO Jv JO I 



Using ( |4.2| ) in (|1.10| ) we have 



1 f°° f°° f Xl f X2 i f Tl i 

Gf(xi,x 2 ;y 1 ,y 2 ) = (-) / dn / dr 2 / Vz\ \ Vz 2 exp—{ dr[{m\ + z 2 ) + 

I JO JO Jyx Jyi Z Jo 

+ jP dr' 2 {m 2 2 + 4 2 )}i(TrPexp^ j^A^z))) = 

j poo roo rx\ rx2 

= (-) 2 / d n dr 2 VzJ Vz 2 expiS qq ~ (4.3) 

2 Jo Jo Jyi Jyi 



S 9q being a kind of effective qq action. Then, in the quenched approximation, keeping 
(lniy)p ert at the lowest order and rewriting appropriately Eq. (|1.4|), one has 

t\nW = i\n(lTrPexpig I dz»A p (z)} = \g 2 H dr[ H dr' 2 D pv {z' x - z' 2 )z^ z' 2 v + 
3 Jv 3 Jo Jo 

+\g 2 f 1 dr[ P d^DM - z'Dz'fzl" + \f f* dr' 2 P dr^D^(z[ - z'[)%z>? 
6 Jo Jo 6 Jo Jo 

+a [ dr[ [ dr 2 5(z[ - z 20 )\z[ - z' 2 | f ds{z' 2 z' 2 - (sz[ T z' 20 + (1 - s)z' 2T z' 10 ) 2 }^ (4.4) 
Jo Jo Jo 



with Zj = dZ3 } T ^ as in Sec. 3. 

8S qq 

Let us now introduce the momenta pj P — — 

j 
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PU = - \ jf dr' 2 D, v {z[ - - t jH dr' 2 D,M - z'[)z'[ v + af^ 



■ i 



| jT dr[D^(z[ - z' 2 )z'r - ~ jT driDM - + af^ (4.5) 



with 



K = sz' 20 z 1T + (1 - s)z[ z' 2T 
f[ = Iz't-z'I / dss s — 

i' 2 = |z' x — z' 2 \ I dss(l — s 



2 



1 u' 



[*3&-u?]* 

fio = - r dr'5{z! XQ - z' w )K - z' 2 | /* ds Z -^ ~ (1 " s) * 2T • U '* 

JO JO 



& = - f - 4o)l< - 41 /* ds^_^k_< (4 . 6) 

1/^10^20 ~~ U s J 2 



Eq.( |4.5|) cannot be inverted in a closed form, however, we can invert it by an expansion in 



2 

ot s — f- and At the lowest order we have 

° 47r m z 



= p f + y dT 2 D^{z[ - 4)p^ + J dr[D^(z[ - z'Dp'Z + aff 
4 M = J# + \g 2 P dr[D, u (z[ - z' 2 )p? + \g 2 P dr' 2 D, v {z' 2 - + aft (4.7) 

O JO O JO 



with 



u 's = SPsoPlT + (1 - s)pi p 2T 

f[ = — \z\ — z' 2 \ ds s 

JO 



1 g: 



/•I u' 

1*2 = — |z x — z 2 | / ds(l — 

JO 



b?o4 - 2 

f / i / _ / I f 1 jP'wP% - (1 - S )P2T • U' s 
J 10 ~~ l Z l Z 2l 



bM-U?]^ 

/j D = K 1 -^|/> ^^^ . (4.8) 

[P10P2O - O 2 







Then we can perform the Legendre trasformation 

f = "Ef^-f^=E P ^[-i(p5» + m 2 )] 

1=1,2 1/0 1=1.2 J ° Z 
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+-g 2 P dr[ P dr' 2 D, v {z[ - 4)pf^ + -g 2 P dr[ P dr?D, u (z[ - z'Dp'fp? 
o Jo Jo o Jo Jo 

+\g 2 P d < P dr 2 D ^2 - 4)&W 

6 Jo Jo 

+a [ dr[ [ dT 2 6(z' w -4o)l z i ~ A\ I Mp'ioPw ~ 
Jo Jo Jo 



(4.9) 



and set 



i=l,2 J0 



(4.10) 



In conclusion at the specified order the qq propagator can be written as 



Gf(x 1 ,x 2 ; yi, y 2 ) = (-) 2 



dn I dr 2 I VziD Pl Vz 2 Vp 2 expi{J2 I dr'^-p'^ + 

Vi 



J-2 
V2 



/ f 2 2 

+-{p j -m j 



4 



g 2 I dr[ I dr^D^z^ - z' 2 )p?p 2 



-g 2 P dr[ P dr'{D,M - z^p? - -g 2 P dr' 2 P dr^D, v (z 2 - 4),%tf 
6 Jo Jo 6 Jo Jo 



3=1,2 



-a / dr[ / dr 2 S(z w - z 2 



2<>J\ Z 1 Z 2 



I' ds{p' 2 p' 2 2 -u 2 '^ + ...}. 
Jo 



(4.11) 



In Eq. ft4.lll) one can proceed as in Eqs. (|2. 12|) — ( p. 17|) and one arrives to the BS equation 
( |3.15| ) in the momentum space with a kernel given by 

d k\d k 2 



I(Pi,P2,p[,P 2 ) = / d z x d z 2 



r 16 2 „ 



■ Z\ + z'y z 2 -\- z' 2 



/IV \ 



(2tt) 8 
)k\k 2 + 



d 4 z[d 4 [ z' 2 e~ i<yPl ~ kl ^ z ' i e ~ i( <p-2-k2)z 



+4a5{ 



Ao + «10 ^20 + 4o m Z l + Z 'l Z 2 + Z 



2 2 71 2 2 

Calculating the integrals in Eq.( f4.12|) , using the definition ( |3.17[) and factorizing the four- 
momentum conservation 5 we arrive to a BS equation of the type ([3.18 ) with a kernel given 

by 0-CD- 

The particular Weyl ordering in Eqs. ( ^4. 12]) , (177), ( |1.8| ) corresponds to have interpreted 
L2f) at the discrete level as 



(4.13) 



Notice that J^ +1 dx^A^x) differs from only i(x^ — Xg^ 1 )A At ( Xa+ ^ 1 ) by terms of the order 
0[(x s -x s „ l ) 3 ] = 0(el). 
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V. SALPETER POTENTIAL AND RELATIVISTIC FLUX TUBE MODEL 



Let us now consider the potential to be used in the Salpeter equation or in the semirel- 
ativistic Schrodinger equation and corresponding to the kernel I given in Eqs. (|1.5|) — (|1.8|) . 

In the scalar model and in the center of mass system, the standard relation occurring 
between the BS kernel and the potential reads 

(k|V|k'> = 1 1 /inst(k,k'), (5.1) 

(2tt) 3 4Jtu 1 (k)iy 2 (k)w 1 (k')w 2 (k') 



where Wj(k) = ymj + k 2 and j inst denotes the so called instantaneous kernel. Precisely J inst 
is obtained from I replacing pj Q and p' by appropriate functions of pj and p'j. The simplest 
prescription corresponds to take 

R0 = p ; = ^ k > + ^ k '> (5.2) 

In the Coulomb gauge the resulting potential is 



,1 It,!, A f 1 4 f i 1 ,, I IK K. I ' LI I , 



2tt 2 3 



1 1 . , , ((k-k')-q 



2 



k-k) 2 q l0 q 20 (k'-ky^ (k'-k) 1 



+ 



1 j rfW *'-^"_L_ [ *» v ^r^ + ^v^lT* 

7T) d 7 2 Om + O20 Cm v (fco v 



(27r) J 7 2 aio + <?2o ?io <?2o 

. 910020 f • i • I^t| m,1 (r o\ 

H — j -(arcsm- + arcsm- -)\ + . . . > (5.3) 

Iqt| |gio| mol J 

with q = and g,o = ^^4r^ ^ and p, = — ; — = . The potential (15.3!) corresponds 

1 1 \J w j (tywj (k) 

to a particular ordering prescription in the hamiltonian ( |1.9|) which, as already noticed, 
is identical to the hamiltonian of the relativistic flux tube model @,[J. Obviously, by an 

2 . 

expansion in ^ one obtains also the semirelativistic potential derived in references 
We notice however that, due to the substitution (5.2) and then to the occurrence of the 
factors P1P2 in (|5.3j), the ordering prescription does not simply coincide with the Weyl 
prescription given in M. 
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VI. CONCLUSIONS 



In conclusion we have established that it is possible to extend the Wilson loop method, 
used to obtain the semirelativistic potential, to the derivation of a Bethe-Salpeter equation. 
We stress that this amounts to show that it is possible to obtain the BS kernel from more 
fundamental arguments than those usually used in the phenomenological application of the 
BS equation. The result is strictly related to the relativistic flux tube model and to the 
semirelativistic potential for heavy quarks. Some additional remarks are in order. 

1. Due to the difficulty in solving directly the BS-equation the general strategy in 
the application of such equation should be this: first solve the three-dimensional 
Salpeter equation for the potential (5.3), i.e. the eigenvalue equation for the 
Hamiltonian (1.9); then evaluate the "retardation correction" by some kind of 



iterative method as it is usually done in the positronium case [ 1 I ] or we have 



done in [12 



2. In Equation (3.12) or (3.17), G2 stands for the complete one particle propagator 
which in principle should be given by (3.5) or its counterpart for the QCD case. 
Due to the absence of a closed path in (3.5), in the present context such quantity 
should be consistently evaluated simply by its perturbative expansion in contrast 
with what sometimes supposed. 

3. The assumption of the approximation (1.4) could seem unjustified in a relativistic 
treatment. Notice however that the important point in QCD is to have some zero 
order extimate of the interesting quantities in a formalism which already provide 
confinement and then to proceed by subsequent corrections. 

4. Finally we notice that the kernel (|1.7]) is highly singular for k' = k, due to the 
occurrence of the factor r under the Fourier transform. This may be inconvenient 
in numerical calculations and makes some equation ill defined. The problem is 
the same that occurres with the linear potential, if one works in the momentum 
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representation, and it is obviously related to confinement. Our philosophy is 
that one should introduce an appropriate infrared regularization (e.g. make the 
substitution r — > re~ Xr = ^^7-) and take the limit A — > only at an advanced 
stage of the calculation. 



APPENDIX A: 

As an example let us derive the contribution in ( |3.14| ) corresponding to the crossed 
diagram in Fig.l. Expanding L\ in (3.11) we have 

Gi(x 1 ,x 2 ;yi,y 2 ) = G 2 (x 1 - y 1 )G 2 (x 2 - y 2 ) + ^v^(^) 2 / dr[ ( dr x [ dr' 2 f dr 2 

A 2 Jo J T [ Jo Jt' 2 

f d 4 z[ f d 4 z 2 r vz x r vz 2 r 1 v Zl r 2 v Z2 d(4 

J J J z' x J z' 2 Jyi Jy2 



exp \{ /J 1 drtt-z? - m? + f dr>{' D{z>{ - z'{')} + [J dr^-tf + J // ^ 



/// 



'1 - " '1 

.2 „ T » 



D{z' 2 > - z' 2 ")}} xexpj{ f 1 dr'( - m\ + | P dr'{' D{z'{ - z'{')] + f 2 ~m 2 2 + 

2 Jo 2 Jo Jo 

+| f dT?D(z> 2 > - 4")] + ^ dr'l f 2 dr»D(4 - 4)} x 
2 Jo 2 Jo Jo 

f 1 + f r ** r 1 - *> + ! r f 



"' - 4") 



1 



+9j ? r dr " ii dr ' w[ -*£)+■••}■ ( ai ) 

The contribution corresponding to the crossed diagram (CD) is that coming from the last 
term in inside the curl bracket in ( |A1[ ). We can write 

x 1 I-X2 rzt rzL pt[ I-T2 

Vzi / Vz 2 / Vz x \ Vz 2 \ dr'i / dr^D{z'{ - 4)D(z[ - z' 2 ) 

[ Jz' 2 Jyx Jy 2 JO Jt' 2 

ml + — dr!" 



exp ^{ jT <[-4' 2 - + I /J 1 - 4")] + j(7 ^<[-4' 2 ~ml + 9 f £ 



2 



D{z>> - 4')]} x exp U f 1 dr'{[-z[ 2 - m\ + f f*' D(4' - 4")] + f * < [-i? - m 2 . + 
2 Jo 2 Jo Jo 



+| jT* - 4")] + / Q Tl dr'l fj dT'>D{z'[ - 4)} 



= -^l) 2 / d4z "J d '4j d 4 z[f d%D{z\ - 4)D(4 - 4) 

roo roo roo roo roo roo 

I dr "l, ** /, dTi l K k K k dT2 

pxi rz[ rz" rx 2 rz' 2 ' rz' 2 

/ Vzi / Vzi / Pzi / Vz 2 / Vz 2 / Vz 2 

J z' x J z' x l Jyi J z'2 J z' 2 J y2 

ex P \{ £ dT>{[-z>? - m\ + I /;" dr1'D{z'[ - 4")] + 

+ £ dr'{'[-z'{' 2 + § \l d<"DW ~ 4")} 
exp \{ £ dr< 2 \-z<? - m\ + | jT" ^"£(4 - 4")] + 
j£ 8 rf<[-if - ml + I /J *f - 4")} x 



'2' " -"2 

rr'/ 



* r f Tl j -iii2 2 1 9l f Tl j mi 7-,/ 111 1111 1 

exp-{ / -m l + — dT x D(z 1 - z x + 

2 Jo 2 Jt{ 

FdTZ'l-Xp - m\ + I d4"D{4 
Jo 2 J T ' 



'" 4" + 



+ ^ r drl' dr? D{z'{' - z%)} x 
2 Jo Jo 

% -{ g l ^ dr'" r" dr""n(7.'" - tS"\ + ^ 



exp 2 if X« dT "' X * dT "" D[z '" ~ z ' n + f ; » d7f / ; ^ fD( 4 " :f 1 



drl' d4'D(4'-4')\. (A2) 



In conclusion, replacing again the last exponential L 2 in ( |A2| ) by 1, we have 

cd = iglgl J d A z[ J d A 4 J d A 4 J d 4 4D(z[ - 4)d{4 - 4)G 2 {x l - 4)g 2 (4 - 4) 

G 2 (x 2 - 4)G 2 (4 - z' 2 )G 4 (zl 4; y u y 2 ). (A3) 

The method can be immediately extended to higher order terms or iterated on (A. 2) ex- 
panding even L 2 . 
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